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>S FOR WIRELESS NETWORKS 



TECHNICAL FIELD 



The field of the invention is wireless communications. The invention is 
particularly concerned with wireless networks, including single hop and multi- 
hop networks. 



Using wireless radio frequency transmissions provides great freedom to 
configure networks. Wireless models become more popular because 
inconvenience and expenses associated with physical connectivity, i.e., through 
wires, are eliminated. The wireless model is also generally more flexible in that 
10 the potential paths for communications are not constrained by any wired 
connection. 

Limitations in the wireless communication model for networks include 
signal interference and power limitations, which both constrain data capacity. 
Signal interference is the natural interference inherent in wireless 
15 communications. Power limitations arise from communication standards, e.g., 
the Federal Communication Commission standards in the United States, that limit 
the power that transmitters may use for a transmission. Power limitations also 
arise from the portable power sources sometimes used by nodes in wireless 
networks. 

20 Power and data rate constraints have been addressed on a node to node 

level, but an efficient and systematic method for accounting for these constraints 
on a network wide basis is believed lacking. Yet, even modest efforts at 
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centralized scheduling in networks have made significant contributions to 
increase network capacity. In the conventional EEEE 802. 1 1 standard there is a 
point coordinate function to centralize scheduling of users in a network. The 
HDR (1EX VDO) standard also schedules transmissions centrally on downlink 
5 channels. There remains a need, however, for a systematic approach that is 
effective and computationally feasible. 



DISCLOSURE OF THE INVENTION 
10 The invention concerns routing, scheduling, and power control methods 

for single and multi-hop wireless networks. A multi-hop network is one in which 

source and destination nodes may communicate directly or through relay nodes. 

Nodes in single hop networks communicate without use of relay nodes. 

Embodiments of the invention may produce an optimal schedule to provide for 
15 the best-case goal for a given parameter. In a preferred embodiment, total power 

is the parameter and total power is minimized for the network. In another 

preferred embodiment, data throughput is the parameter, and throughput is 

maximized for the network. 



20 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1. is an exemplary network topology for application of methods of 
the invention and used in a numerical simulation of an embodiment of the 
invention; 

FIG. 2 is a graph of total average transmit power versus network 
25 throughput for three cluster schedules X, Y and Z reached in a numerical 
simulation and applied to the network topology of FIG. 1; 

FIG. 3 is a graph of total average transmit power versus network 
throughput for the cluster schedules X, Y and Z reached in a numerical 
simulation and applied to the network topology of FIG. 1 ; 
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FIG. 4 is an example transport style network for application of methods of 
the invention and used in a numerical simulation of an embodiment of the 
invention; 

FIG. 5 plots the ratio of throughputs achieved by scheduling in accordance 
5 with the invention compared to TDMA with increasing values of ambient noise 
for the network topology of FIG. 4; 

FIG. 6 plots the total average power with increasing traffic rates for the 
session between node 1 and node 5 in the network topology of FIG. 4; 

FIG. 7 illustrates an example asymmetric diamond topology for 
10 application of methods of the invention and used in a numerical simulation of an 
embodiment of the invention; 

FIG. 8 plots throughput ratios for simulations on the FIG. 7 network 
topology; 

FIG. 9 plots average power required to achieve higher traffic demands for 
15 the example FIG. 7 network topology; 

FIG. 10 illustrates an example hierarchical network topology for 
application of methods of the invention and used in a numerical simulation of an 
embodiment of the invention; and 

FIG. 1 1 plots the ratio of the total maximum throughput for a routing 
20 policy produced by application of methods of the invention versus that using only 
minimum energy paths as a function of ambient noise and in the example 
topology of FIG. 10. 



25 BEST MODE OF CARRYING OUT THE INVENTION 

The invention concerns routing, scheduling, and power control methods 
for single and multi-hop wireless networks. A multi-hop network is one in which 
source and destination nodes may communicate directly or through relay nodes. 
Nodes in single hop networks communicate without use of relay nodes. 

30 Embodiments of the invention may produce an optimal schedule to provide for 
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the best-case goal for a given parameter. Iq a preferred embodiment, total power 
is the parameter and total power is minimized for the network. In another 
preferred embodiment, data throughput is the parameter, and throughput is 
maximized for the network. 

5 Embodiments of the invention also provide an integrated routing, link 

scheduling and power allocation policy for a general multi-hop network that 
minimizes the total average power consumption to support minimum average rate 
requirements per link. Multi-hop networks implementing methods of the 
invention can support higher throughputs than with conventional approaches to 

10 radio resource allocation. Methods of the invention may be implemented in 
multi-hop networks that are time synchronized between transmitters, and are 
characterized by channel conditions that remain essentially constant over several 
time slots. 

Methods of the invention are generally applicable to single and multi-hop 

15 networks, but will not perform as efficiently in networks with rapidly changing 
channel conditions. The invention is particularly well-suited for application to 
networks of wireless nodes in a small geographical area, such as a WI-FI 
network, or an IEEE 802.11 style network. 

According to another aspect of the invention, large-scale networks are also 

20 handled. For large-scale networks, a hierarchical approach is applied whereby 
the set of links is partitioned into clusters. The links in each cluster are scheduled 
somewhat independently from the links in other clusters, and intercluster 
interference is modeled as static ambient noise. Interaction between clusters is 
modeled with a fixed-point equation, which determines the level of inter-cluster 

25 interference. At a higher level, scheduling can be done at the cluster level, rather 
than at the link level. 

For embodiments that minimize total power, we find our optimal power 
control policies are such that each node is either transmitting at peak power to a 
single receiver or not transmitting at all, which makes it simpler to implement 

30 from the standpoint of computational resources and speed. This also means that 
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nodes can go to sleep much of the time, when they are not transmitting or 
receiving, potentially conserving a great deal of energy. The optimal link 
schedule time-shares a small number of optimal subsets of links (~ L + 1 for a 
network with L links) in order to achieve the required data rates. The collection 
5 of optimal transmission modes depends on the level of ambient noise and the 
capacity allocations on links. 

A byproduct of solving the scheduling and power control problem yields 
link sensitivities that provide the change in the nrinimal total average power with 
respect to a perturbation of the required data rate on a link. The minimal required 

10 average power is a convex function of the required data rates on each link. This 
leads us to a routing algorithm. 

Testing of the invention revealed that that using minimum-energy paths to 
route traffic and TDMA (time division multiple access) to schedule links is near- 
optimal in low noise regimes and in the case where the required average data 

15 rates are small. As the level of ambient noise increases, or the required average 
data rates increase, our optimal routing policy exploits the use of multiple paths 
to support higher data rates, and the optimal scheduling policy makes use of 
concurrent transmissions on links. We find that non-minimum energy paths can 
be exploited to increase throughput, despite the fact that all links share a common 

20 bandwidth. 

Preferred embodiments of the invention will now be described. Example 
implementations and simulations are provided, but do not limit the invention in 
its broader aspects and are not necessary for understanding or implementing the 
methods of the invention as fully and separately described. Methods of the 
25 invention may be implemented in any node of a network, or in a dedicated 
centralized control node. 

Suppose there are N stationary nodes, labeled by the integers 1,2, . . ., N. A 
set e of L e = | 8 1 transmission links, among the possible N (N - 1) links between 
nodes, constitutes a network topology. 
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For simplicity of exposition, we divide time into slots, each of equal 
duration and indexed by the positive integers. Transmissions begin and end on 
slot boundaries. For a given link I = (i,j), the transmitter node i uses a signal 
power P m (l) in slot m. The path gain from node i to node j is given by G(i,j), and 

5 can model the effects of signal attenuation due to distance, channel fading and 
shadowing, as well as antenna gain patterns. We model path gains G(i,j) as being 
constant over the time frame being scheduled. As mentioned, the invention is 
most effective with relatively constant channel conditions. If the channel 
conditions are rapidly changing, the path may still be assumed constant, but the 

10 time period of scheduling might be kept smaller, e.g., scheduling over 100 time 
slots instead of over 200 in response to changing channel conditions. 

The transmitting and receiving nodes of link I are denoted by T(l) and R(l) 
respectively. The received signal power at node R(l) in slot m from the 
transmitter T(l) is thus given by P m (l)G(T(l), R(l)). However, signals emanating 

15 from other transmitters appear to the receiver R(l) as interference, and there is 
ambient thermal noise to contend with as well. The signal to interference and 
noise ratio (SIR) for link / in slot m is defined as 



20 



G{T(l) 9 R(l))P m (l) 

7mW = 



Ejbdb* P m {k)G(T(k), R(l)) + n m (1) 



where nj is the "ambient noise" power at the receiver at node j. For now 
25 we assume that the ambient noise power at each node does not change from slot 
to slot, but we will relax this assumption later in a hierarchical approach that 
represents a preferred embodiment. 

The data rate of the single wireless channel embodied by link / during slot 
m can be modeled as a function of the SIR, y m (I). It is assumed that the achieved 
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data rate on link 1 during slot m, X m (l), is a linear function of the SIR, y m (I), and 
the total available bandwidth W. For simplicity of explanation, all links are 
assumed to share the same frequency band of width W. Using the Gaussian 
approximation to compute the bit-error-rate (BER), the data rate of a link / with a 
5 tolerable BER of 10"* and using BPSK modulation is given by 



Let Pm = (Pm(l),Pm(2),...Pm(Le)) be the network power vector for slot m. 
Let P" 1 "*^ ) be the maximum transmission power for node L Let s(0 be the links 
15 in s that originate at node L 

Each node must conform to the peak transmission power constraint in 
every slot: 



X m {l) = W Jm (l) where W = 



W 



2g-log 2 10" 



(2) 



10 



0< ^Pm(Z)<P n,ax (0and 




Using (2) and (1), the achieved data rate for link I in slot m is 



20 



Xm=W'(; 



G(W),R(l))P m (.l) 



St : it * iG(T(k), R(l))P m (k) + mm 



The long-term average rate of link I is then defined as 



(Z) = lim-fp&(0. 



For each link I, let C(l) be a given minimum acceptable average data rate, 



i.e. 



we must have 
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(4) X avs (I) > C(Z), for all I e e . 

Define the required minimum average rate vector as 
C = (C(l),C(2),...,C(Le)). 

The average power consumed by the transmitter for link Z, P™ 8 ^), is 

l m 

P™9(1)= Urn — T^flbCO- 



Define the average network power vector as 

p«* = pan p«* (2)j m ^ pavg (Le)) 

There may or may not exist a sequence of network power vectors Px,P2,... 
that satisfy (3) and (4). If there does exist a sequence of such network power 

— -avg 

10 vectors, our aim is to ininimize the total average power given by h(P ), where 

(5) h(p ms ) = Y, paVS W- 

lee 

A systematic approach is used by the invention to compute a 
schedule of network power vectors Pi t Pi f ... , satisfying the peak power 

-—avg . 

constraints (3) and the minimum average rate constraints (4), such that h(P ) is 
15 minimized. This primal problem is succinctly stated as 

(6) min{fc(F v * )} subject to (3) and (4). 

Note this optimization involves choosing optimal power levels in 
each slot for each transmitter. Using a duality approach the problem is reduced to 
a convex optimization problem over a single slot. The duality approach is 
20 preferred to solve the stated minimization in view of constraints. More generally, 
however, a linear program may be used to reach a solution, as will be appreciated 
by artisans. 

Duality Approach 

Let n = (%<!), »r (2) , -,n R(U) ) be ambient noise power vector. Let the value of 

25 the 
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optimal cost in the problem (6) as a function of C and n be denoted by 
H(C,n) . If for a given value of C and n , no schedule of network power vectors 
exists satisfying (3) and (4), define H(C,n) =+<~. It can easily be seen that for a 
fixed value of C, H(jcji) is monotone increasing in n . Also, it can be shown that 
5 for a fixed value of n , H{C, n) is a convex function of C . 

Define a set of dual variables 0 = {0<X),02,...,0(L k )} for each link. Define 
a "potential" function V according to 



Note the absence of the minimum rate constraints (4) in the definition of 
g(fi). A primary advantage of the restatement is that the number of variables to 
15 be optimized is smaller. A secondary advantage is that the optimal values of the 
dual variables yield link rate sensitivities, used in embodiments of the invention 



V(P,C, X) = KP) +£/?(Z)[C(Z) - X(0] • 



tee 



10 



The dual objective function is defined as 




for routing. For any non-negative vector p , using (4), note that 



H (C, n) :> min {v<P"" , C, x"* )} , subject to (3) and (4) 




(7) 



= g(#>. 



20 



Thus, for any non-negative vector J3 , the minimum value of the objective 
function in the optimization problem (6) is lower bounded by g(p) . This 
observation leads to the dual optimization problem, 

,„x max{<?(£) :0>O}. 
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A convexity argument can be used to show that in fact we have 




Computation of g(/3) involves optimizing over all schedules of 



network power vectors satisfying the peak power constraint in every slot. 
5 However, since the potential function v(P,C,X) is linear in P,C, and X, it 
follows that g(j3) can be computed by an optimization over a single slot, i.e. 

q(P) = min{vCP,C, X) : Pe S P ] , 
where X = X(P) is defined as in (2) and (1), and S P is defined as the 
polyhedral set 

10 S P = {p : Z tee(0 P(Z) < P™* (0 for all i and 0 < P(l) for all I e e} . Let M be the 

number of extreme points of S P , and let Kf , m = 1, 2, ...M denote the extreme 
points of S P . By definition, each point in S P can be represented as a convex 
combination of Kf , m = 1, 2, .. M . An upper bound on M is 2 U . If Mi is the 
number of links emanating from node i, then M = II^Cl+M,) ■ 
15 A convexity argument then shows that 



Thus, the complexity of computing g(j3) is O(M). In fact, note that 
V(Pm ,C,X(Pm)) is an affine function of /?, which implies that the dual 
objective function g(/8) is the pointwise minimum of M affine functions, i.e. 



ascent algorithm converges in a finite number of iterations. Alternately, by 
introducing a slack variable, the dual optimization problem of equation (8) can be 





(10) 



20 



An iterative ascent algorithm may be used to solve the dual optimization 
problem of equation (8). Exploiting the piecewise linearity of g(0) , the iterative 
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reduced to a linear program and be efficiently solved using interior point 
methods. 

We'll now discuss a termination condition of the iterative ascent 
algorithm. For a given n , if the required minimum rate vector C is infeasible, 
5 i.e. H{C,n) = +oo , it can be seen that g(fi) is unbounded from above over the set 
of all /? > 0 . Otherwise, it is easy to see that 

Thus we can safely declare C and n infeasible if we find a value of J3 
10 such that 

Computing an Optimal Link Scheduling Policy 

Solving the dual problem (8) yields the optimal dual variable vector and 
extremal network power vectors P* such that 

The vectors P* ,£ are extreme points of S P that represent optimal 
transmission modes and the power levels of the transmitters in those transmission 
modes. In other words, assuming H(C,n) is finite, an optimal schedule of 
network power vectors exists that consists solely of these extremal network 
power vectors P*'. This means that in the optimal policies we find, in each slot, 
every node is either not transmitting at all, or is transmitting at the max i mum 
possible peak power to exactly one receiver. Moreover, since a hyperplane in 
(4 +1) dimensional Euclidean space is determined by 4 + 1 linearly independent 
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points contained within it, an optimal policy can be constructed that consists of at 
most 4+1 extremal network power vectors P*'\ Let Kbe the number of 
extremal network power vectors P*' 1 such that (1 1) holds. Define 
K = min {K, L e + 1} . We thus have 

5 

for all i such that l<i< K . 

In any optimal policy (assuming H(C,n) is finite), the average rate on 
each link I, is exactly equal to C(l). Otherwise, we could decrease the 

total average power further without violating any rate constraints, thus 
10 contradicting the optimality of the solution. Let X*/(l) denote the rate of link I 
corresponding the optimal extremal network power vector P*\ Assuming 
H(C,n) is finite, it can be shown that (9) guarantees the existence of a "weight 

vector" X ={A*(X),A\2),A\k )}such that the equations (12) and (13) below are 

satisfied: 

15 (12) 

^A*(i)X^(J) = C(l) for all links le£ 

5=1. 

ic 

V>*0") = 1 where A*(s) > 0 . 

^ (13) 

Note that K < L £ +1 , so the above system of equations is always overdetermined. 
There must be K +1 of these equations which are linearly independent. Once we 
20 have identified such a set of equations, we solve them by inverting a 
(K + l)x(K+l) matrix. 

The value of X*(i) indicates the relative frequency at which the extremal 
network power vector P* 1 is utilized in an optimal policy. This then determines a 
family of optimal policies, which can be constructed by alternating between the 
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extremal network power vectors in accordance with the weight vector X . A 
specific optimal policy, for example, can be constructed using a pseudo-random 
number generator to select which extremal network power vector is used in each 
slot, such that the network power vector P* 1 is used with probability 

5 

Reducing Complexity 

The algorithm developed to solve the dual problem considers all M 
transmission modes. However, in the optimal schedule we have no more than 
L £ +1 transmission modes. Since M can be as large as 2 Ls , the complexity of 

10 minimizing the dual objective function can be exponential. Thus, reducing the 
set of possible transmission modes considered can greatly reduce complexity. A 
natural consequence of the invention for a network of nodes with omni- 
directional antennas is that no node transmits or receives data at the same time in 
the optimal schedule due to high self interference. This observation can reduce 

15 the number of candidate transmission modes significantly. In general, in many 
cases it is possible using heuristics to eliminate many obviously inefficient 
transmission nodes from consideration, since they are unlikely to be used in an 
optimal policy. 

To summarize the preferred scheduling embodiment just discussed and 
20 make some general observations about it, artisans will appreciate that the general 
steps included are to measure the channel parameters between transmitter and 
receivers in order to assess the corresponding data rates on each link for each 
possible transmission mode. Using a primal optimization approach, duty cycles 
of each transmission modes in an optimal scheduling policy can be calculated by 
25 formulation of the problem as a linear program. Alternatively, using a dual 
optimatization approach, we can first identify a relatively small set of 
transmission modes that comprise an optimal scheduling policy. The 
identification of this small set of transmission modes can also be posed as a linear 
program. Once such a set of transmission modes is identified, a set of linear 
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equations can be solved to identify the duty cycles of each transmission mode in 
the set. Once an optimal scheduling policy is calculated, it can be disseminated to 
all transmitters and receivers so that the optimal transmission modes can be 
realized in a synchronous manner. 

5 

Hierarchical Scheduling 

For scheduling a network with a large number of links, increased speed 
may be realized by taking a hierarchical approach while applying the general 
scheduling discussed above. This can avoid a large amount of computation when 
10 the number of links renders the scheduling calculation too expensive. The basic 
idea is to schedule links hierarchically.. In the hierarchical approach, we 
partition the set of edges s into / disjoint subsets of edges s(l), e(2),... e(J), 
called clusters. We thus have e = Uj^e(j) . Since the links in the same cluster 

must coordinate to implement the scheduling algorithm discussed above, it is 
15 convenient to group links into clusters such that the links belonging to the same 

cluster are in geographic proximity to each other. 

Let C (j) be the minimum rate vector for cluster j, i.e. the components of 

C (j) consist of the components of C which are the minimum rates for each of the 

links belonging to cluster j. Similarly, let n (j) be the ambient noise vector for 
20 clustery, i.e. the components of n (j) consist of the components of n which are 

the ambient noise powers for each of the receivers belonging to clustery. 

Decoupling the Scheduling for Each Cluster 

When we partition the links in the entire network into clusters, we assume 
25 that the number of links in any given cluster is limited so that if the network 

consisted of only the links in the cluster, it is computationally feasible to use the 
method described above to schedule the links in the cluster. Let Hj(C(j)Mj)) 
be the minimal total average transmit power required for meeting the average 
data rate requirements for only those links belonging to clustery. In other words, 
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the links in all other clusters besides cluster j are assumed to be absent or idle in 
the definition of Hj(C(j),nU)) . The value of Hj(C(j)Mj)), as well as a 
corresponding schedule for the links in cluster j 9 can be found by application of 
the embodiments described above. 
5 Of course, the noise at a given receiver will not only consist of 

interference from other links in the same cluster and ambient thermal noise, but 
interference from links in different clusters will also be present. The power of the 
interference from other clusters will vary from slot to slot, in accordance with the 
link schedules used in each of the clusters. However, suppose we model the 

10 additional noise from interference in other clusters as static ambient noise. 

Specifically, suppose that N(j) = (^(l),iV y (2),...,iV y (L W) )) consists of the 
total "average ambient noise" for the receivers of links in cluster j. In order to 
calculate the value of N(j), first define ^(ik) = (iT , (lXir , ( 2 ).--^(^))) as 
the average transmission power vector for the transmitters of the links in cluster 

15 k. The value of P™* (k) is obtained as a byproduct of computing H k (C(k),n(k)) 
and the associated schedule for cluster k. Indeed we have 

H*(c(*).^*»=z^jr(0- 

Consider link V in cluster t The transmitter of link /' transmits with 
average power P™ 8 (/') . Thus the average power of the interference caused by 
20 link V in cluster k to the receiver of link I in clustery is G(T(l') 9 RQ))P^ 8 (/') . The 
value of Nj(l) is obtained by summing the total contribution from all such links V 

in all other clusters other than j, and adding the thermal ambient noise at receiver 
n R(l) , i.e. we have 

(14) 
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Let N be the concatenation of the noise vectors N(\),N(2),...,N(J). From 
the discussion above, the vector N is determined by the global rate requirement 
vector C and the global ambient thermal noise vector n . We can denote this 
mapping by a function F(.,.), i.e. 
5 (15) N = F(C,n). 

Coordinating the Clusters: Finding a Fixed Point 
In general we will have ~N > n , so we need to take into account the effect 
of inter-cluster interference in the scheduling of each individual cluster. The 

actual ambient noise vector we seek is N , where N satisfies the fixed point 
10 equation 

(16) N* = F(C,N*). 

This equation may be solved iteratively. In particular, let be the 
estimate of the solution at the tf h iteration. We can start with a noise vector that 
consists only of ambient thermal noise, i.e. N° =n. We then find a sequence of 

15 estimates using the recursion = F(C,N k ). We have performed a series of 
numerical experiments, and have found in every case that such an iterative 
algorithm finds a fixed point. 

If each component of F(C,N) is a monotone function of N , it is easy to 
establish that the iterations above will also converge to the smallest solution of 

20 equation (16), assuming a solution exists. In the case where each cluster is 
sufficiently distant from all other clusters, all transmissions emanating from a 
given cluster can be modeled as a point source. In this case F(C,N) is a 
monotone function of H j(C(f),N(j)) for all./, which in turn are monotone 

functions of N . Hence, the convergence of the iterative algorithm can be proved 
25 when all clusters are sufficiently distant from each other. We also believe, but 
cannot prove, that F(c,N) is a monotone function of N in all cases. 
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To summarize, we propose that each cluster first calculates a schedule, 
assuming all other links in other clusters are idle and there is no inter-cluster 
interference. After this calculation, we can compute the incident average 
interference at each receiver from links that belong to other clusters. Taking into 

5 account this average inter-cluster interference, we again calculate a schedule for 
each cluster. Proceeding inductively, we recalculate the incident average 
interference for each receiver from the transmitters belonging to other clusters. 
We then again calculate a schedule for each cluster using the new value of 
incident average inter-cluster interference, and so on. 

10 Accounting for Non-constant Inter-cluster Interference 

After a fixed point solution N* to (16) is found, we have a transmission 
schedule for each cluster in which we have assumed that the power of inter- 
cluster interference is constant The constant is the average inter-cluster 
interference obtained from the solution to the fixed point equation (16). 

15 However, since the emitted power from a cluster varies from slot to slot, 

this assumption is not true. Fortunately, essentially due to the convexity of the 
function (a+ bx)" 1 with respect to x, the actual average SIR for any given link 
will be at least as large as the average SIR obtained by assuming that the power 
of inter-cluster interference is constant. Thus, the actual achieved data rate for 

20 each link is at least as large as it would be assuming constant power for inter- 
cluster interference. In practice, each cluster could communicate its schedule to 
nearby clusters, and the data rate used in each slot can be appropriately adjusted 
to account for the non-constant inter-cluster interference. 

To prove this, suppose link I belongs to clustery and let X"* (I) be the 

25 actual average data rate achieved for link I, taking into account that the power of 
the inter-interference varies from slot to slot. Also, let X*™ 8 (I) be the average 
data rate achieved link I assuming that there is constant power ambient noise 
N*(l) . We wish to show that X°" (!) > X^ (I) . 
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Recall that N*j(l) is determined from the solution to the fixed point 
equation (16). In particular, if n R(i) (m) is the total average ambient noise incident 
at the receiver R(l) in slot m, then 

1 n 

Nm= lim -Y^n m {m) . 
(17) ra=l 
Suppose the fraction of time in which the transmission mode i is 
used in the schedule for cluster j is where transmission mode i is given by 

Pj' = (P/''(1),C( 2 ). • ^ S »W be set of slots with index less 1:11311 
or equal to n where transmission mode i is used in clustery. In fact we assume a 
stronger condition than in (17), namely 



(18) 

for each transmission mode i in the schedule for clustery. This will be 
true, for example, if the pseudorandom number generators for selecting 
transmission modes in each of the different clusters are independent. Moving on 
to the proof, we have 



X™*{1) = lim -yzx m {i) 

77S=L 
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where the first inequality follows from (18) and the convexity of the 
function (a+fot)" 1 with respect to x. 

Scheduling at the Cluster Level 

5 Above, we have assumed that the links in each cluster can be active in 

each slot, i.e., all of the clusters are active in each slot. In certain situations, this 
may result in very high levels of inter-cluster interference, particularly when the 
links in different clusters are in close proximity. 

As an alternative, we can group the clusters into different edge sets. A 

10 given edge set consists of the union of some subset of all the clusters, and each 

cluster belongs to exactly one edge set. The set of all slots is partitioned in such a 
way that each slot is reserved for links that belong to a given edge set. 

With a suitable grouping of clusters into edge sets, we can insure that the 
clusters belonging to a given edge set are sufficiently separated in space. This 

15 will then limit the inter-cluster interference for clusters belonging to the same 

edge set. Each edge set can then be scheduled using the method discussed in the 
previous subsections. 

More generally, we could formulate the cluster scheduling problem in a 
similar way to the link scheduling problem considered a above on an entire 

20 network basis. However, it appears that satisfactory schedules can be found by 
only considering a small collection of edge sets, where each edge set is assigned 
the same fraction of the overall system bandwidth. Moreover, by limiting the 
number of edge sets to a small constant, the scalability of the overall approach is 
significantly enhanced, both in terms of computation complexity, and the ease to 

25 which the overall algorithm can be implemented in a distributed way. We also 
note that, alternatively, the edge sets can be multiplexed in frequency rather than 
in time, as artisans will recognize. 

Routing an d Capacity Allocation 
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As described above, scheduling conducted according to the invention can 
find an optimal link scheduling and power control policy to support a given 
traffic rate on each link in the network while consuming minimum power. 
Additionally, with a change in the statement of the constraints, scheduling of the 
invention may be used to find optimal link and power control policies to 
maximize data capacity. We also have developed a method to construct a policy 
which allocates traffic rates on each link in accordance with a routing algorithm, 
given a traffic matrix which describes the rate at which traffic is to be carried 
between specific source-destination pairs. The routing algorithm is guided by the 
computation of the minimum cost associated with optimal scheduling and power 
control for supporting traffic rates on each link that have already been allocated, 
as described in the previous section. 

For example, consider a series of iterations, where at a given iteration we 
wish to support a small amount of additional traffic from one of the given source- 
destination pairs. Initially, no traffic is supported, i.e. the required minimum data 
rate on each link is zero. We compute the optimal values for the dual variables, 
g(j3) Using g(/?) = minjv(P, C, X) : Pe S P ] , we can estimate the cost of 
supporting additional traffic on each link using a first order Taylor series 
approximation. For example, suppose we wish to evaluate a given route r. Let r 
be the L E dimensional vector, such that the component r(l) = 1 if link / is 

included in the route r , and r(l) = 0 otherwise. If we route an additional e units 
of traffic along route r, the required data rate vector will increase from 

CtoC+er. The cost of supporting the additional e units of traffic along route r 
can be determined from the approximation 

H (C =e r) = H(C)+g 20* (I)KD • 

tee 

Therefore, in order to determine the best possible path along which an 
additional increment of traffic should be routed, we can use a shortest path 
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algorithm (such as Dijkstra's algorithm or the Bellman-Ford algorithm), where 
we set the weight for link / to the sensitivity value 0* (I) . After we have 
identified the shortest path, we attempt to support the additional traffic along this 
path. If the new value of C is infeasible, we either reduce the step-size e until 
5 the additional traffic can be supported, try supporting additional traffic for 
another source-destination pair, or otherwise reach a prespecified termination 
condition. 

Indeed, since #(•) is convex, we can use the optimal routing framework 
developed by Bertsekas et al. (e.g., see, D. Bertsekas and J. Tsitsiklis, "Parallel 

10 and Distributed Computation: Numerical Methods", Prentice Hall 1989.). In this 
class of algorithms, any feasible routing is first chosen. By incrementally 
adjusting the routing in an iterative manner, flows are shifted towards shortest 
paths at each iteration, and eventually a globally optimum routing is obtained. 
There are many possible variations, for example regarding whether flow is 

15 shifted simultaneously for every source destination pair, or whether flow is 

shifted one at a time for each source destination pair. In the present context, we 
are faced with a more difficult problem, since #(•) is not differentiate 
everywhere. In addition, since evaluating the function is not trivial, it is 
more difficult to ascertain whether a given routing is feasible, i.e. the capacity of 

20 a given link depends on the rate carried by other links. Nevertheless, our initial 
numerical experiments and the framework developed by Bertsekas et. al. suggest 
that such an approach may be quite useful in practice. 

The potential energy efficiency of a given link can be evaluated, assuming 
that all other links are idle. For our model, the minimum transmission energy 

25 required per bit on a link (i,j) is inversely proportional to G( ij). A minimum 

energy path is defined as a path which has the smallest possible energy required 
to transport a fixed amount of data, irrespective of the delay associated with 
transporting that data. Previous work on routing in wireless networks (for 
example see, S. Singh, M. Woo and C. S. Raghavendra, "Power-aware routing in 
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mobile ad hoc networks", Proc. of Fourth Annual ACMREEE Conference on 
Mobile Computing and Networking pp. 181-190, Oct. 1998) is focused on energy 
efficiency without regard to delay, i.e., such work proposes using minimum 
energy paths. 

We find that the paths over which traffic is allocated in some optimum 
routings does not always correspond to using minimum energy paths. Our 
simulations revealed, however, that when required data rates are sufficiently low, 
our routing framework will apparently always use minimum energy paths 

Experimental Simulation Examples 

To illustrate the invention, we will now discuss some simulation examples 
in the next section. 

Hierarchical Scheduling (Example 1^ 

We consider a large fully connected wireless network consisting of 80 
nodes and 150 links as shown in FIG. 1, distributed over a 200m x 200m field. 
The network was divided into 16 clusters labeled by integers 1, 2,... 16, and each 
cluster consists of 5 or 6 nodes and 8 or 9 links. A link belongs to a cluster if the 
transmitting node is within or on the boundary of that cluster. Moreover, each 
link in the network belongs to exactly one cluster. The minimum average rate on 
all links in the network is set to the same value. The peak power of all links is 
identically set to 1 Watt and W = 10 6 . We applied our hierarchical approach over 
the entire network by dividing it into edge sets in 3 different ways denoted by 
cluster schedules X, Y and Z. Cluster schedule X schedules a single edge set 
concurring all 16 clusters. Cluster schedule Y schedules two edge sets each 
consisting of clusters in a checker board pattern, i.e., { 1,3,6,8,9,1 1,14,16} and 
{2,4,5,7,10,12,13,15}. Each edge set in Y is scheduled for half the time. Cluster 
schedule Z schedules four edge sets containing clusters { 1,3,9,11 }, {2,4,10,12}, 
{5,7,13,15} and {6,8,14,16} for a fourth of the time each. We plot the total 
average transmit power consumed by all the links in the network as a function of 
the (minimum) network throughput for each cluster schedule in FIGs. 2 and 3 
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corresponding to different values of thermal noise n = 0.05m W and n = 0.2 mW 
respectively. The thermal noise at all receivers is the set to be the same. 

We find that for very low data rates (between 0 to 41 Mbps in fig. 3) it is 
energy efficient to enable all the 16 clusters simultaneously, regardless of the 

5 value of thermal noise. As the data rate along links increases (between 41 to 78 
Mbps in FIG. 3), cluster schedule Y is the most energy efficient of the three 
configurations considered. For data rates beyond 78 Mbps (in FIG. 3), the energy 
efficient strategy is to enable the cluster schedule Z. For both values of thermal 
noise, Z is capable of supporting higher throughput than that achievable by X and 

10 Y. 

For low rates, the optimal strategy in each cluster is a simple TDMA of 
links with long idle times. Therefore, the intercluster interference is small even 
when clusters are next to one another, making cluster schedule X the most energy 
efficient. However as the data rate requirement on links increases, links in each 

15 cluster are active for longer durations resulting in higher intercluster interference. 
In this case, schedules Y or Z are more energy efficient. We observed similar 
results to that in FIG. 6 for randomly generated network topologies as well. 

Observe that by selecting a schedule based on the rate requirement we can 
realize an energy versus rate graph which is the minimum of the curves 

20 corresponding to the cluster schedules X, Y and Y in FIG. 6. In fact, by suitable 
timesharing of the cluster schedules we can achieve any rate in FIG.6 which lies 
in the convex hull of the minimum of the three cluster schedule curves. 

Our iterative approach converges quickly to the fixed point solution for all 
cluster schedules considered making it practical for implementation in large 

25 wireless networks. We also observed that as the distance between the clusters in 
an edge set increases, the difference between the actual data rate and that 
obtained by assuming constant inter- cluster interference (over a slot) decreases. 
This is attributed to the weaker dependence of the transmission strategy of each 
cluster to the intercluster interference. For cluster schedule Z, we found that the 
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actual data rate was within 8% of the value obtained using constant inter cluster 
interference when the thermal noise n = 0.2 mW. 

Example 2 String Topology 

5 

A string topology consists of a row of nodes connected by means of 
directed links as shown in FIG. 4. This topology is useful as a transport network 
to carry data through short hops but over potentially long total distances through 
relay nodes. Node 1 is the source of data and node 5 is the sink. Each node in this 
10 example has an omnidirectional antenna, and all nodes have identical peak 
transmission power constraints of 1 Watt. The path-loss between nodes is 

modeled as an inverse square of the distance, i.e., G(i, j) = * . . The ambient 

noise power at all the nodes is assumed to be constant. We compare the 
maximum throughput of the single session between node 1 to node 5 achieved by 

15 our link scheduling policy to a TDMA scheduling policy. The TDMA policy 

schedules exactly one link at a time, activating each link at maximum power for a 
fourth of the time. We plot the ratio of throughputs achieved by the our link 
scheduling algorithm to TDMA with increasing values of ambient noise in FIG. 
5. The ambient noise is normalized with respect to the peak received power (from 

20 a single transmitter) in the network and is scaled by log 2 (.) (for this and the 

remaining examples). The maximum throughput achieved by our policy is higher 
than TDMA for all values of ambient noise. In low noise regimes, TDMA is a 
near-optimal policy. However, as the level of ambient noise at nodes increases, 
the gains in throughput of our policy increase significantly. The optimal policy 

25 schedules multiple transmissions at the same time, even in fairly low ambient 
noise regimes, thereby exploiting spatial reuse and thus outperforming TDMA. 
This effect is more pronounced in moderate and high ambient noise regimes 
where the optimal policy schedules concurrent transmissions even though they 
are in close geographic proximity. 
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To get a sense of how energy efficient the optimal policy is as it supports 
increasing traffic demands, we plot the total average power with increasing traffic 
rates for the session between node 1 and node 5 in FIG. 6. For this experiment, 
we fix the level of the logarithm of the normalized value of ambient noise to - 
5 0.67. This value corresponds to the 3rd data point on the x-axis in FIG. 5. We 
also set W= 10 7 . 

FIG. 6 indicates that the total minimum average power is a linear function 
of the traffic demand for average rates below 4 Mbits/sec. Below 4 Mbits/sec, our 
scheduling method functioned essentially as TDMA. Specifically, only one link 

10 is active at any time, each link is active for the same fraction and for some 

fraction of time all links are idle. As the traffic demand of the session increases 
beyond 4 Mbits/sec, the total average power increases at a faster rate. For traffic 
load between 4 Mbits/sec and 4.98 Mbits/sec, the transmission modes in the 
optimal policy are given by [{(1,2)}, {(2,3)}, {(3,4)} and {(1, 2), (4, 5)}]. For 

15 traffic loads above 4.98 Mbits/sec to the maximum throughput possible, the 
optimal transmission modes are [{(1,2)}, {(2,3)}, {(3,4)}, {(1,2),(3,4)} and 
{(1,2),(4,5)}]. We see that in order to support high traffic loads, the optimal link 
scheduling policy activates a large number of links simultaneously. Note that the 
transmission modes in the optimal policy are half-duplex. Therefore, 

20 transmissions in the optimal schedule cannot consume a total average power of 
more than 2000 mW for this topology. 

Example 3 Diamond Topology 

We also consider an asymmetric diamond topology as shown in FIG. 7 to 
25 illustrate the substantial increase in throughput by splitting traffic over multiple 
routes including paths that are energy inefficient. Nodes in this topology are 
equipped with omnidirectional antennas, and the peak trans- mission power each 
node is fixed at 1 Watt. Node 1 is the only source of data and node 4 is the sink. 

We assume G(l,4) = — — T ; the path loss between all other nodes is given by the 
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inverse square law of distance. We consider our integrated routing and 
scheduling algorithm over this diamond topology, where all possible routes from 
the source to the destination are allowed except for the single hop route (1,4). For 
comparison purposes, we also considered our scheduling algorithm where only 
the links (1,2), and (2,4) are allowed to be used, which corresponds to routing 
over a single nrinimum energy path 1 — > 2 — ► 4. 

We compare the ratio of throughputs for these polices in FIG. 8. Clearly, 
using multiple paths yields higher throughputs for all values of ambient noise. 
The increase in throughput by splitting traffic over multiple paths is significant 
even for moderate levels of ambient noise. This result is somewhat surprising 
since all links are sharing a common bandwidth and we have a per-node peak 
power constraint. An intuitive explanation of this is that by splitting the traffic 
over both paths, the transmission modes {(1,2),(3,4)} and {(2,4),(1,3)}, which do 
not have much interference can be alternated in time. This results in substantial 
delivery of data to the destination in every slot. If we only allow the path 1 — ► 2 
-> 4, in order to avoid self interference, the links (1,2) and (2,4) must be active 
one at a time, so that node 4 will only be delivered data in every other slot. 

To illustrate energy efficiency issues, we plot the total minimum average 
power with increasing traffic demand at node 1 in FIG. 9. The logarithm of the 
normalized value of ambient noise for this point is 0.45 and corresponds to the 
3rd data point on the x-axis in FIG. 8. We set W = 10 7 . The maximum data rate 
achieved by TDMA using the minimum energy path 1 — ► 2 — > 4 is 3.65 Mbit/sec. 
In comparison, our policy is capable of supporting a data rate of up to 5.05 
Mbits/sec by using both the minimum energy path as well as the non- min i m um 
energy path. For traffic loads below 3.65 Mbits/sec, our method produces a 
policy that functions as TDMA over the path 1 — » 2 — > 4. As the traffic load 
increases beyond 3.65 Mbits/sec, the method results in an optimal policy that 
starts using both paths by scheduling transmission modes {(1, 2), (3,4)} and 
{(1,3), (2,4)} in addition to TDMA transmission modes. The average power 
required to achieve higher traffic demands beyond 3.65 Mbits/sec increases with 
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a greater slope as can be seen in FIG. 9. The maximum throughput of 5.05 
Mbits/sec is achieved by splitting traffic by using both paths and scheduling 
transmission modes {(1,2), (3,4)} and {(1,3), (2,4)} for a dominant fraction of 
time and scheduling transmission modes {(1,2)}, {(2,4)} for the remaining time. 
5 The optimal policy transmits a fair amount of traffic on both available paths, but 
transmits a greater share of its traffic on its minimum energy path. We have noted 
a 10 to 20% reduction in the throughput gains of the optimal policy than what is 

shown in FIG. 8 for the case when G(l, 4) = — ^— - . 

d(lA) 

Example 4 Hierarchical Topology 

10 A hierarchical topology as shown in FIG, 10 represents a complex 

topology such as an access network, where the number of nodes in each tier 
decrease as we get closer to the access point node 6. This experiment highlights 
the significant gains in throughput by jointly routing and scheduling to an 
alternative policy that optimally schedules links over pre-determined minimum 

15 energy routes. Nodes 1, 2 and 3 are sources of data and node 6 is the sink. We 

assume G(l,6) = — ^—t and G(2,6) = — T ; the path loss between all other 

d(l,6) d(2,6) 

nodes is given by the 

inverse square law of distance. The minimum energy paths for the source 

nodes are: {1 3 6}, {2 — ► 3 — ► 6} and {3 — ► 6} respectively. Our integrated 

20 routing and scheduling algorithm is allowed to utilize all possible paths to route 

their data to node 6. We constrain the rates of each source to be identical. We plot 

the ratio of the total maximum throughput of our routing policy to that achieved 

by using only the minimum energy paths, as a function of ambient noise in FIG. 

11. 

25 We find that at low noise levels, our policy supports a throughput at least 

30% more than the data rate achieved using minimum energy paths alone. As the 
level of ambient noise increases, our policy supports a significantly higher 
throughput. The base policy schedules links in a TDMA sequence for all values 
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of ambient noise. The optimal policy schedules concurrent links even for low 
levels of ambient noise. In the optimal policy, nodes 1 and 2 route their traffic 
over non- minimum energy paths 1 — > 4 — > 6 and 2 — > 5 — > 6 respectively, apart 
from the minimum energy paths, for all values of ambient noise. For low ambient 
5 noise, the optimal policy schedules the transmission modes {(1,3)}, {(2,3)} and 
{(3, 6)} according to TDMA but also allows modes {(1,4),(5,6)} and {(2,5),(4,6)} 
to be active for reasonable fractions of time. As the level of ambient noise 
increases, the optimal policy splits node 3*s traffic over paths 3 — > 4 — > 6 and 3 
— > 5 —> 6. In fact, for sufficiently high ambient noise, the optimal policy 
10 schedules a large number of concurrent links, the dominant transmission modes 
being {(1,4), (5,6), (2,3)} and {(2,5),(4,6),(1,3)} and {(1,4), (2,5), (3,6)}. 

Generalized Embodiments 

Methods of the invention as illustrated in the specific embodiments are a 

15 teaching of more general embodiments of the invention, as artisans will 

appreciate. As one example, methods of the invention can minimize the total 
average transmitter and receiver energy required to achieve the target capacity on 
links in the network. The constraints and the value being minimized/maximized 
is simply changed. Another constraint we could impose is a maximum peak 

20 power constraint for an arbitrary collection of nodes, rather than per each node. 
This may be useful to limit signal emissions from an entire network or 
subnetwork. 

Minimizing Total Average Transmitter and Receiver Power 
25 The total energy required to transmit a packet over a slot is the energy 

expended in the electronic circuitry to transmit at the desired power level. The 
efficiency of power amplifiers is typically around 10% and therefore, the total 
power consumption at node T(l) link is roughly 10P(l). Energy involved in 
coding (e.g. convolution coding) packets may be treated as negligible. 
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The energy required to receive a packet of information is comprised of 
two components. The first component is the electronic energy required to operate 
PLL's, registers, etc. and this is proportional to the data rate of the link. The 
second component is the energy required to decode the packet. While the 
decoding energy increases greatly as the SINR of the link decreases, this 
component of energy can be substantially lowered (by orders of magnitude) by 
using Application Specific Integrated Circuits to decode the received signal. 

We therefore ignore the decoding energy component and coarsely model 
the power expended by a receiver R(l) as a linear function of the data rate of link 
Z. Let P R (l) denote the power expended by node R(l) to receive data at rate X(l). 
Then P R (l) = Q(l)X(l) where 0(7) is the energy expended per bit that is constant for 
each receiver node R(l). The new cost function that we would like to minimize is 

h(p aV8 )= 2>« pcv * (0+ZW(0 d9) 

lee lee 

15 where h(P ) is the sum total of average transmission and receiver power. 

The dual objective function corresponding to the new primal function (19) would 
have an additional average receiver power term. The dual problem corresponding 
to the new primal problem is similar to (8) except that without loss of generality, 
we can constrain f3(l) > 0(1) for each link Ze s. This constraint, 0(1) > 0(1) , allows 

20 us to exploit the concavity property and use the link scheduling and power 

control algorithms originally developed to solve problem (8), for the new dual 
problem. 

Accounting for receiver power in our method changes the choice of routes 
and the rate allocations on links in each route. The changes in the routing policy 
25 are more pronounced in a network where nodes are in fair proximity to each 
other. Using multiple hops in such topologies would consume higher receiver 
energy than using a single hop, making multi-hop routing inefficient. In such 
topologies, ignoring receiver energy consumption in the problem formulation 
would have yielded multi-hop routing as the optimal routing policy. 
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As described above, embodiments of the invention provide an integrated 
routing, scheduling and power control framework that is well suited for slow 
fading wireless channels which are relatively constant for long durations of time. 
Outdoor wireless channels with strong line-of-sight paths typically tend to have a 
5 very stable channel strength and could be approximated as a constant for a time 
interval (block fading model) ranging from a few to several minutes. By 
periodically measuring the channel and accounting for the changing traffic 
requirements for each source-destination pair, the invention can compute the 
optimal routing, scheduling and power control method and adapt it to the 
1 0 changing network conditions . 

While specific embodiments of the present invention have been shown and 
described, it should be understood that other modifications, substitutions and 
alternatives are apparent to one of ordinary skill in the art. Such modifications, 
substitutions and alternatives can be made without departing from the spirit and 
15 scope of the invention, which should be determined from the appended claims. 

Various features of the invention are set forth in the appended 

claims. 



